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NUMMARY 

An addition to the computer program of NASA SP-273 Is given that permits 
transport property calculations for the gaseous phase. Approximate mixture 
formulas are used to obtain viscosity and frozen thermal conductivity. Reac- 
tion thermal conductivity Is obtained by the same method as In NASA TN 0-7056. 

Transport properties for 154 gaseous species have been selected for use with 
the program. 


INTRODUCTION 

An option for calculating the dilute-gas transport properties of complex 
chemical mixtures has been added to the chemical equilibrium program of 
reference 1. Computer programs for calculating mixture transport properties 
are presently available (e.g., ref. 2). Reference 2 gives transport properties 
in the form of collision Integrals and uses the Hlrschfeider ''rigorous" method 
(ref. 3) for the viscosity and thermal conductivity of nonreacting mixtures. 

The Brokaw method (refs. 4 and 5) Is used for the contribution of chemical 

thermal, conductivity. The -objectives are (1) to use a method for 
calculating nonreacting (frozen) mixture properties which Is simpler than that 
presented In reference 2, While stlU retaining the same accuracy, and (2) to 
use a simpler representation for the properties of pure species and their 
J]" ary 1n J erac tlons . These objectives shoul substantially reduce calculation 
times or Increase the size of the system which can-be calculated economically. 

new program transport properties for pure species and several binary 
Interactions are represented In the functional form suggested In reference 6. 
(See the section "Transport-Property Dataset.") Mixture viscosity Is Obtained 
by an approximate method In a form first proposed by Sutherland (ref. 7). Mix- 

form f^r^oncl! 1 / s °5 ta1ned By an approximate method In the 

form first proposed by Wasslljewa (ref. 8), which Is s-lmllar to Sutherland^ 

form for viscosity. Reaction thermal conductivity Is obtained In the same 

d V?k' re * erence 2 * lhe total thermal conductivity of the mixture Is 
the sum of the frozen and reaction thermal conductivities. 

*h 1h1 *!” epor t does not repeat the material presented In reference 1 (such as 
the equations and numerical techniques for obtaining chemical equilibrium com- 
positions, mixture thermodynamic properties, and various dDDllcatlons of fhA«« 
proles), only *MU1ona1 ■tSrlll reUtl^to 


APPROXIMATE MIXTURE METHODS 
Sources of Approximate Methods 

In addition to the "rigorous" method for obtaining the mixture transport 
properties mentioned In the INTRODUCTION, numerous approximate methods may be 
found In the llteratu Some approximate methods for mixture viscosity are 
described. In references 7 and 9 to 14, while some approximate methods for 
mixture frozen thermal conductivity may be found In references 8 , 12 , and 15 
to 20. Most of these approximate methods have forms similar to those first 
proposed by Sutherland for viscosity (rof. 7) and by Wassl ljewa.for thermal 
conductivity (ref. 8 ), which can be written as 
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where 


number of gaseuus species for transport calculations 


mole fraction of species 1 

viscosity of species 1 

n m 1 x viscosity of mixture 

thermal conductivity of species 1 

X m 1 x,fr frozen ~ thermal conductivity of mixture 

viscosity Interaction coefficient between species 1 and J (<p * <p ) 

♦u frozen- thermal conductivity Interaction coefficient between species 1 
1 and J ( 4 ^ * 4 »jj) 


The approximate methods of this type differ only In the expressions 
H> 1 j and 4^ j • For example, Hlrschfelder *s first approximation to his 

"rigorous" method for mixture viscosity contains an expression for the 
In equation ( 1 ) which Is used by some other authors as a starting point 
deriving other approximations to q>,,. Hlrschfelder 's expression (ref. 
eq. (8.2-26)) Is equivalent to ” 
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where Is the molecular weight of species 1, Is a ratio of collision 
Integrals (defined In ref. 3, eq. (8.2-15.)), and Is a quantity which can 

can be defined In several ways, some of which are dlscus-sed later. Various 
expressions for derived by other Investigators from equation (3) or by 

Independent analysis can be put In the following general form for the pur- 
pose of comparison: 

*1J * k 1j ro 1J (4 
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and where r Is an exponent and k^j Is a function having several forms. 

Two examples are given to Illustrate equation (4). For the first example, 
Brokaw (ref. 12) gives an expression for «p,. which Is equivalent to the 
following: J 
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In this case r * 1 and 
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As a second example, Wilke (ref. 11) gives the following expression for q>^: 
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and r » 1/2. 



A summary of for several approximate methods Is given In table I. 

These methods differ from one another only In the exponent r.. The quantity 
which appears In equation (6) can be obtained In several ways as follows, 

An expression for tu, Is given by Hlrschfelder (ref. 3, eq. (8.2-21)): 
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where 8^’ Is a collision Integral and a function of c^, and and 

are Lennard-Jones Interaction-potential parameters for species 1 and J. 
For 1 » j, equation (10a) reduces to the following form for pure species: 
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where V - Q i • Hlrschfelder suggests 

the following approximations for a;.d c^: 
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An alternative expression for can be derived based on the following 

approximation for 8^., given by Svehla and McBride (ref. 2, eq. (36)): 
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The expression for derived from equations (10a), (10b), and (13) Is 
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where 


Is defined by equat1on.(9) . Substitution of n. 


eqyat1on_(15) Into equation (6) gives the exact expression for <p^ In 

equation (8). Thus, If the approximation given In equation (13) Is used, 
equations (6) and (8) are Identical. 

In addition to equations (10a) and (15), values of n,, may also be derived 
from experimental data on mixture viscosity. These three J methods for obtaining 
generally give different numerical values. Therefore, even the "rigorous" 

method gives different values for mixture viscosity depending on what value Is 
used for 

Expressions for (refs. 12 and 15 to 19) are given In table II. The 

expressions from references 17 and 19 are the same as that from reference 18 
except for the factors of 1.065 and 0.866. The formulas of Lindsay and Bromley 
(ref. 15) contain parameters known as Sutherland constants which can be esti- 
mated In various ways. Sutherland constants for this report were obtained from 
viscosity by means of the following formula (derived In ref. 21); 
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The Sutherland Interaction constant S^ was obtained from the following 
approximation: 




(16b) 


In addition to the foregoing methods there Is an "averaging" method for 
thermal conductivity due to Burgoyne and Weinberg (ref. 20). A generalized 
version of this method Is 
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An analogous formula for viscosity Is 
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Selection of Approximate Method- 


A number of mixture methods were tested to discover the extent each repro- 
duced experimental viscosity and thermal conductivity data for several binary 
systems. The methods tested were the "rigorous" method; the approximate 
methods, whose expressions for <p^ and appear In tables I and II; and 

the averaging method of reference 20. All the methods reproduced the experi- 
mental data by better than 3 percent for at least some of the test cases. 
However, each method also produced errors exceeding 7 percent In at least one 
case. 

An example of the spread In accuracies ts shown In table III for the Ar-NH 3 
system using the experimental data of reference 22. The experimental data are 
for eight different compositions of Ar and NH 3 at a temperature of 308 K. Cal- 
culations were made with values of from equations (10a) and (15) and a 

value of = 175 derlyed from experimental data. In order to derive a~ 

value from experimental data, It Is necessary to select the form of the mixture 
equatloa with which It Is to be used. For this example, the mixture method of 
reference 12 was selected with k^j from equation (7) and with the exponent 
r * 1. As expected table III shows that the best results vaverage absolute 
error = 0.2 percent) were obtained using the method of Brokaw (ref. 12) and 
the derived from experimental data with the same method. 

For the binary systems tested, comparisons with experimental data showed 
that no one method of calculation could be judged superior to all others In all 
cases. The selection of a particular method represents a compromise between 
overall accuracy and ease of computations. To obtain viscosity of mixtures, 
the expression for from equation (6) Is used In equation (1) If values 

of are available In the dataset discussed In the section "Transport- 

Property Dataset". Otherwise from equation (8) Is used (which Implies 

using a value of from eq. (15)). 

Brokaw's mixture method for frozen thermal conductivity (ref. 16) was 
selected Inasmuch as It reproduced the experimental data for the test cases 
slightly better than the other methods tried. The same values used for 

for viscosity are also used to calculate for frozen thermal conductivity. 

Reaction thermal conductivity Is obtained by the same method as used In 
refereftce 2. 


COMPUTER PROGRAM FOR TRANSPORT PROPERTIES 

The computer program was written In FORTRAN IV and tested on an IBM 
370/3033 Computer at the Lewis Research Center. A logical variable TRNSP1, was 
added to NAMELIST/INPT2/ of reference 1. Transport calculations are carried 
out only If TRNSPT = T has been Included In the namelist Input. Table IV lists 
the I/O (Input/Output) units used by the program. A transport-property data- 
set, five subroutines, and BLOCK DATA are Involved In transport-property 
calculations. 

One of the five subroutines, SEARCH, has two purposes. One purpose Is to 
search the unformatted thermodynamic data on I/O unit 4 for all possible 
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species In any specified chemical system and to store thermodynamic coeffi- 
cients for these possible species In a COMMON block. A second purpose, If 
TRNSP1 » T, Is to search the unformatted transport data on I/O unit 8 for all 
possible species In any specif led. chemical system and to store transport coef- 
ficients for these species on I/O unit 3. Subroutines TRANP and INPUT and the 
entry 0UT4 In subroutine 0UT1 are used exclusively for transport-property cal- 
culations. BLOCK DATA contains headings for transport property output as well 
as other Information required by the program. The fifth subroutine, UTRAN, Is 
used only for converting formatted transport properties to an unformatted form. 


1 ransport- Property Dataset 

Transport properties for 154 pure species and several binary viscosity 
Interaction quantities were least- squares fitted to the following form, 

used Inreference 6: 


= A In 1 + ~ ~ + D (19) 

1 r 


The coefficients were generated to give viscosity In units of pP and thermal 
conductivity In units of pW/cm K. The temperature range was divided Into two 
Intervals to be consistent with the same Intervals selected for thermodynamic 
properties, namely, 30G to 1000 K and 1000 to 5000 K (although the program 
permits any three Intervals for any species). Therefore, In the present data- 
set, each pure species or binary Interaction may have four sets of coeffi- 
cients: two sets for viscosity (high- and low-temperature Intervals) and two 

sets for thermal conductivity (high- and low-temperature Intervals). If only 
viscosity or only thermal conductivity Is present, only the two sets of coeffi- 
cients for that property are given. The code word TRAN must precede the 
formatted sets of coefficients, and the code word LAST must follow. The format 
for the dataset containing these coefficients Is given In table V. 


In n 
In x 


Subroutine UTRAN 

Coefficients for generating thermodynamic and transport properties are used 
by the program In an unformatted form In order to reduce computer time consid- 
erably. Conversion of the coefficients for transport properties from formatted 
to unformatted form Is done by subroutine UTRAN. The formatted coefficients 
are read In as Input on 1/0 unit 5. Subroutine UTRAN converts the formatted 
coefficients and writes them In the unformatted form on I/O unit 8. After 
conversion, the coefficients are available for future runs In unformatted form 
on I/O unit 8, and subroutine UTRAN need not be used again until the formatted 
dataset Is revised. 


Subroutine INPU1 

After composition has been determined for a specified thermodynamic data 
point (such as an assigned temperature and pressure), subroutine INPUT 
selects coefficients for a maximum of the 50 most abundant gaseous species for 
that point from among the transport coefficients previously stored on 1/0 
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unit 3. Depending on the accuracy desired, the number of species used by the 
program may often be reduced by means of an optional input parameter TRPACC. 
If, for example,. TRPACC Is set equal to 0.98, the program will consider only 
those species which add up to 0.98 of the total molar composition (adding In 
the order of largest to smallest species mole fraction). The default value of 
TRPACC Is 0.99995. The. mole fractions, of these selected species are first 
normalized to give a total of 1.0 before mixture transport properties are cal- 
culated. The transport properties are calculated from the selected coeffi- 
cients at the current temperature. 


Subroutine TRANP 

Subroutine TRANP sets up the equations and solves for mixture viscosity, 
frozen thermal conductivity, equilibrium thermal conductivity, and equilibrium 
specific heat. The method and equations for obtaining equilibrium specific 
heat and the reaction contribution to thermal conductivity are given In 
reference 2. 


Entry 0UT4 In Subroutine 0UT1 

Entry 0UT4 handles the output of the transport-property calculations. A 
variable format Is used which Is the same as the format used for thermodynamic - 
propertles output (described In ref. 1). 


DESCRIPTION OF PROGRAM OUTPUT 


A sample problem containing transport-property output Is given In table VI. 
Transport properties are given for mixtures of gases only. Therefore, for mix- 
tures containing condensed species, two values of equilibrium specific heat are 
given. One value (which Includes gaseous and condensed-species contributions) 
Is given with the output labled THERMODYNAMIC PROPERTIES, while a second value 
(which contains only gaseous-specleS contributions) Is given with the output 
labled TRANSPORT PROPERTIES. For each point a maximum of the 50 most abundant 
gaseous species Is used for calculating mixture transport properties. Thermal 
conductivities are given for the assumptions of both equilibrium and frozen 
gaseous compositions. Viscosity Is always given In units of mP. An option Is 
provided In the namelist Input for units of thermodynamic properties. Units 
for thermal conductivity are cons-lstent with the option selected for thermo- 
dynamic units. That Is, thermal conductivity Is In units of mW/cm sec If 
SlUNIT * T Is Included In the Input dataset and mcal/cm K sec otherwise. The 
Prandtl number has the usual definition of Cpn/x, while the Lewis number as 
defined In reference 2, Is 


Veac c p.fr 

x fr c p,reac 


( 20 ) 
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SELECTION OF TRANSPORT PROPERTIES 


Sources of Properties for Pure Species 


Transport properties for all species except those discussed In the follow- 
ing paragraphs were taken from the compilations of Svehla (ref. 23).. Tor 
several species thermal conductivities were calculated for this report from 
viscosities by means of the following equations, given In reference 2: 


31 R_ 

tr,1 * 4 M 1 n 1 


( 21 ) 


X 1nt,1 s M 1 n 1 
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11 C 1nt,1 " 



X 1 s X tr,1 + X 1nt,1 


( 22 ) 

(23) 


where 


. § 

1nt,1 ‘ R " 2 


(24) 


and 


^D. rot. 1 
rot,1 = R 


(25) 


For atomic species, ^ - 0 and = x^ r 4. C ro t \ Is assumed to be 1 
for linear species and l!5 for nonlinear species. 


Argon. - References 6, 24, and 25 contain critical reviews of viscosity 
data for Ar. References 6 and 25 include more recently available experimental 
viscosity data In their evaluation than found In reference 24. Additionally, 
recommended values from references 6 and 25 are In close agreement (less than 
1 percent difference from 200 to 2000 K). By contrast, the values of 
reference 24 agree with values from reference 25 at room temperature but differ 
by 8.5 percent at 2000 K. viscosities and thermal conductivities were taken 
from reference 25 in order to have a consistent set of values. 


Atomic carbon . - Collision Integrals are given In reference 26 for every 
1000 K from 1000 to 25 000 K. These collision integrals were extrapolated 
down to 300 K and converted to viscosities. Thermal conductivities were 
obtained by means of eauatlon (21). 

Me thane . - Viscosities selected by references 6, 24, and 27 to 30 are all 
In extremely close agreement (less than 1 percent difference). The tempera- 
ture range for viscosity In these references Is from 200 to 1000 K, although 
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not every reference covers the entire range. The average of viscosities In 
these six references was selected. Thermal conductivities were selected to be 
the average of the values from references 27 and 30 to 32. 

Carbon monoxide . - Viscosities from references 24, 27, 33, and 34 were 
reviewed. The values selected by reference 27 are about 2 to 3 percent lower 
than those of the other three references In the range from 500 to 1000 K and 
were not used. The averaged values from the other three references are 
extremely close to the N 2 viscosity values of reference 35. The reference 35 
N 2 values were therefore selected as tiie viscosity for CO Inasmuch as the N 2 
values are given to 2000 K. Thermal conductivities were selected to be the 
average of the values In the four references 27, 31, 33, and 36. 

Carbon dioxide . - Viscosities from the five references 6, 24, 27, 37, and 
38 were reviewed. Only those of references 6 and 38 were used to obtain aver- 
age values. The values In reference 37 were not used because, according to 
figure 3 In that reference, the calculated values appear to be about 0.5 to 
2 percent too high from 400 to 800 K. The older values of references 24 and 
27 were not used because they appear to be about 4 percent lower than the more 
recent values In references 6 and 38. Thermal conductivities were averaged 
from those given In the five references 27, 31, and 39 to 41. 

Atomic hydrogen . - High-temperature viscosities (above 1000 K) were cal- 
culated from the collision Integrals given In reference 18. Low- temperature 
viscosities (200 to 500 K) were taken from reference 42. Values from 500 to 
1000 K were Interpolated. Thermal conductivities were calculated by means of 
equation (21). 

Molecular hydrogen . - Viscosities of H 2 were taken to be the average of the 
values In the four references 6, 24, 27, and 43. Thermal conductivities from 
the four references 27, 31, 40, and 43 were reviewed. The selected values are 
as follows: from 200 to 600 K, average values from the four references; from 

1100 to 5000 K, reference 43 values; from 700 to 1000 K, values "faired" 
between the other selected values. 

Water. - Viscosities aid thermal conductivities from 373 to 1073 K were 
taken from reference 44 and extrapolated to higher temperatures by means of 
equations given In reference 45. 

Atomic nitrogen . - The high-temperature (1000 to 15 000 K) collision 
Integrals given In reference 46 were converted to viscosities and extrapolated 
to lower temperatures. Thermal conductivities were obtained by means of 
equation (21). 

Ammonia . - Viscosities were selected to be the averages of those given In 
references 23, 24, and 47. Thermal conductivities were selected to be the 
average of those In references 27, 31, 48, and 49. 

Nitric oxide . - The two sets of viscosities and thermal conductivities for 
NO from references 2 and 24 are fairly close. The reference 2 viscosities were 
selected Inasmuch as a least-squares fit of these values gave less error than 
a similar fit of the reference 24 values. The reference 2 thermal conductiv- 
ities were selected In order to be consistent with the viscosities. The 
selected transport properties are also given In reference 50. 





Nitrogen dioxide . - A problem In obtaining transport properties for NO2 lies 
In the fact that there Is considerable association of NO2 to N2O4 In the tem- 
perature region of about 300 to 400 K. Thus experimental data are open to 
Interpretation as to whether measurements are f.or the pure species, an equili- 
brium mixture of NO2 and N2O4, or some other combination of the two species. 
References 24 and 51 Interpret the experimental data of reference 52 differ- 
ently. There Is close agreement for temperatures above 400 K (where the equi- 
librium mixture Is 98 percent NO2) but considerable difference otherwise. 
Inasmuch as reference 51 obtained viscosity data for NO2 and N2O4 simultane- 
ously from the- reference 52 values, those results were selected for viscosity. 
Thermal conductivities were obtained from viscosities by means of equations 
(21) to (23) with Z rot * 

Di atomic nitrogen . - A comparison was made among the critically evaluated 
and selected viscosities given In the three references 6, 24, and 35. The 
agreement between the data from references 6 and 35 Is excel leat (only a few 
tenths of a percent difference in the entire temperature range from 200 to 
2000 K). Reference 24 Is in good agreement with the- others at. low tempera- 
tures but differs by as much as 8 percent at 2000 K. Viscosities and thermal 
conductivities were taken from reference 35 Inasmuch as reference 6 does not 
contain thermal conductivities. 

Nitrous oxide . - Viscosities from references 2 and 24 for N2O are very 
close (less than 0.6 percent difference from 200 to 1000 K). The reference 2 
values were selected since they are given to 5000 K. The thermal conductiv- 
ities from reference 31 appear to be closer to the experimental data than do 
those In reference 2 and also closer to one recent value given In reference 53. 
Ihe reference 31 thermal conductivities were selected up to 700 K. Above 700 K 
they were calculated from viscosities by means of equations (21) to (23) with 
Zrot a *• 

Atomic oxygen . - Two sets of high-temperature viscosities were compared 
(refs. 27 and 46). These agree to within 2.5 percent from 3000 to 5000 K. 

The viscosities of reference 46 were selected since these are given from 1000 
to 5000 K, whereas the reference 27 values are given only above 3000 K. Low- 
temperature viscosities (200 to 300 K) were taken from reference 54. Viscosi- 
ties from 400 to 900 K were Interpolated. Thermal conductivities were obtained 
by means of equation (21). 

Molecular oxygen . - Viscosities from references 6, 24, 27, and 35 were 
reviewed. The averages of these viscosities are within 1 percent of the 
reference 35 values. The reference 35 values were selected. Thermal conduc- 
tivities from references 27, 31, 35, 55, and 56 were reviewed.. The averages 
of these agree to within 1 percent with reference 35 from 200 to 600 K, but 
differ more at higher temperatures (3 percent from 1000 to 1300 K). According 
to Hanley (ref. 35), "...most of the data for the thermal conductivity coeffi- 
cient for both nitrogen and oxygen seem unreliable outside the range of 150 to 
600 K . " Thermal conductl vltes were also taken from reference 35 to be consis- 
tent with the selected viscosities. 

Hydroxyl radical . - fcxpei lifiental transport data for the OH radical are not 
available. The collision Integrals estimated In reference 57 and the rota- 
tional collision number Zfot * 8 estimated In reference 2 were used to 
calculate viscosities and thermal conductivities. 
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Sources of Properties for Binary Interactions 

The transport-property dataset currently provides coefficients for obtain- 
ing for only three pairs. These coefficients were taken directly from 

reference 58. Values of for other pairs are estimated In the program by 

means of equation (15) using k^ from equation (9). 


APPENDIX - SYMBOLS 


A, B, C, D 


* * 


A 11» A 1j 


C 1nt,1 

C rot,1 


c 

D 

D 

k 


P 

1 

U 

U 


Le 


M 

m 


1 

U 


n 


R 


r 



1 


coefficients In eq. (19) 

ratio of collision Integrals, dimensionless 

constant-pressure heat capacity, J/kmol K 

term defined by eq. (24), dimensionless 

term defined by eq. (25), dimensionless 

constant-pressure specific heat, J/kg K 

self-diffusion coefficient, m 2 /sec 

binary diffusion coefficient, m 2 /sec 

function defined by eq, (7) or (9), dimensionless 

Lewis number, defined by eq. (20), dimensionless 

molecular weight of species 1, kg/kmol 

molecular weight function defined by eq. (5), dimensionless 

number of gaseous species Included In transport calculations, 

dimensionless 

universal gas constant, 8314.41 J/kmol K 
exponent In eq. (4) 

Sutherland constant for specie* 1, K 

Sutherland Interaction constant for species 1 and j, K 

temperature, K 


*1 

Z rot 

e 1 

C U 

n 

n i 

n u 

V 

K 


U 

d 1 

v 1j 

♦1J 


«(! 


?,2)' 


mole fraction of species 1, dimensionless 

rotational collision number 

depth of potential energy well for species 1,1 

quantity defined by eq. (12) 

viscosity, kg/m sec 

viscosity of species 1, kg/m sec 

quantity defined by eqs. (10a) and (15), kg/m sec 

thermal conductivity, W/m K 

thermal conductivity for species 1, W/m K 

quantity In table II, w/m k 

molecular diameter for species 1, A 

quantity defined In eq. (11), A 

viscosity Interaction coefficient (eq. (6) or (8)), dimensionless 

!lJ?* n ;ft ,er !?? 1 ‘ co ? duct1vUy interaction coefficient (eq. (2) and 
table II), dimensionless ' 

collision Integral, dimensionless 
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ft 1j 

collision cross section defined In eq. (14), *2 

Subscripts: 

fr 

frozen contribution 

1 J 

Index for species 

Int 

Internal contribution 

mix 

for the mixture 

reac 

chemical reaction contribution 

rot 

rotational contribution 

tr 

translational-contribution 
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TABLE IV. - INPUT/OUTPUT DATA FILES 


I/O 

unit 

Data file 

3 

Scratch unit 

(1) Used during conversion of formatted to unformatted data 

(2) Unformatted coefficients to generate transport properties for a 

particular chemical system 

4 

Unformatted coefficients to generate thermodynamic properties for all 
species 

5 

Formatted input. Case input (described in ref. 1) and formatted 
coefficients to generate thermodynamic and transport properties: 
(Note: Formatted coefficients need to be processed only once. 

Program will write the unformatted files 4 and 8.) 

6 

Output 

8 

Unformatted coefficients to generate transport properties for all 
species and pairs of species 
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TABLE VI. - THERMODYNAMIC EQUILIBRIUM PROPERTIES AT ASSIGNED TEMPERATURE AND PRESSURE 



